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ABSTRACT 

We present the results of two extensive Rossi X-ray Timing Explorer observations 
of large X-ray flaring episodes from the high-mass X-ray binary pulsar LMC X-4. Light 
curves during the flaring episodes comprise bright peaks embedded in relatively fainter 
regions, with complex patterns of recurrence and clustering of flares. We identify pre- 
cursors preceding the flaring activity. Pulse profiles during the flares appear to be simple 
sinusoids, and pulsed fractions are proportional to the flare intensities. We fit Gaussian 
functions to flare peaks to estimate the mean full-width-half-maximum to be ~68 s. 
Significant rapid aperiodic variability exists up to a few hertz during the flares, which 
is related to the appearance of narrow, spiky peaks in the light curves. While spectral 
fits and softness ratios show overall spectral softening as the flare intensity increases, 
the narrow, spiky peaks do not follow this trend. The mean fluence of the flare peaks 
is (3.1 lb 2.9) X 10^*^ ergs in the 2.5-25 keV energy range, with its maximum at ~1.9 x 
]^q41 gj-gg_ The flare peak luminosity reaches up to (2.1 it 0.2) x 10^^ ergs s~^, far above 
the Eddington luminosity of a neutron star. We discuss possible origins of the flares, 
and we also propose that inhomogeneous accretion columns onto the neutron star polar 
caps are responsible for the observed properties. 

Subject headings: accretion, accretion disks — pulsars: individual (LMC X-4) — stars: 
neutron — X-rays: bursts — X-rays: stars 



1. INTRODUCTION 

Out of ~200 known neutron star /black hole X-ray binaries, more than 40 sources have been 
observed with a wide range of coherent spin period of ~0. 07-1400 s, identifying the central source 
as a rotating, highly magnetized neutron star powered by accretion (van Paradijs 1995; Bildsten 
et al. 1997). These accretion-powered X-ray binary pulsars (hereafter. X-ray pulsars) often show 
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abrupt increases in X-ray luminosity known as "bursts" or "flares." About 30 % of the X-ray 
pulsars have been observed to show bursts and/or flares, including (1) the low-mass X-ray binary 
4U 1626-67 (McClintock et al. 1980), (2) the bursting pulsar GRO J1744-28 (Fishman et al. 1995; 
Kouveliotou et al. 1996), (3) the high- mass X-ray binaries such as LMC X-4 (e.g., Levine et al. 
1991; Woo, Clark, & Levine 1995; Levine, Rappaport, & Zojcheski 2000; Moon & Eikenberry 2001) 
and SMC X-1 (Angelini, Stella, & White 1991; Moon, Eikenberry, & Wasserman 2002), and (4) 
the transient Be binary system EXO 2030+375 (Parmar et al. 1989). Although they arc easily 
identified as sudden increases of X-ray luminosity in most cases, the detailed properties of the 
bursts/flares differ much from source to source. For instance, whereas bursts from the high- mass 
system SMC X-1 have been detected only once despite extensive studies of the source ^, another 
high-mass system LMC X-4 has been observed with large, semi-regular flaring episodes over the 
last two decades by various X-ray satellites. This diversity in bursts/flares suggests that many 
different mechanisms are responsible for the observed phenomena. 

Most theoretical studies of the origin of bursts/flares from X-ray pulsars have relied upon one 
of three following mechanisms: (1) accretion disk instabilities, (2) nuclear burning on the surface of 
the neutron star, and (3) non-uniform stellar winds from mass-donating companions. Obviously the 
most well known case for the accretion disk instability is the X-ray bursts from the bursting pulsar 
GRO 1744-28, for which a thermal-viscous instability in the accretion disk has been investigated to 
explain the origin of the bursts (e.g., Cannizzo 1996, 1997). The X-ray bursts from GRO 1744-28 
share many properties with the Type II X-ray bursts from the Rapid Burster (Lewin et al. 1996) 
and the X-ray flares from SMC X-1 (Moon et al. 2002). On the other hand, the presence of a 
He-rich dwarf companion led Brown &: Bildsten (1998) to suggest the possibility that --^lOOO-s flares 
from 4U 1626-67 are caused by carbon burning under the neutron star surface. Taam, Fryxell, &: 
Brown (1988) attributed the origin of the flares from EXO 2030+375 to an accretion disk instability 
caused by the inhomogeneity in the high-velocity (~550 km s""-*^) winds from a Be star companion. 

LMC X-4 is a persistent, disk-fed, high-mass system with a pulsational and an orbital period 
of ~13.5 s and ~1.4 d, respectively (Kelley et al. 1983). It is one of a few X-ray binaries showing 
a super-orbital period (~30.3 day; Ilovaisky et al. 1984) possibly caused by a precessing, tilted 
accretion disk. The optical companion is a 14th magnitude 0-type star of ~15 Mq (Sanduleak 
& Philip 1977). Over the last two decades, LMC X-4 has exhibited large X-ray flaring episodes 
which have made it one of the most regular and representative flaring X-ray sources. During the 
flares, the X-ray spectrum softens considerably, and the pulse profiles become simple, sinusoidal 
(Levine et al. 1991, 2000). Recently, Moon &; Eikenberry (2001) reported two relatively long-time 
scale structures in the flares: quasi-periodic oscillations of ~0. 65-1. 35 and ~2-20 mHz. They also 
found that the amplitudes of the flares are almost exactly modulated by those of the neutron star's 
coherent pulsations (or vice versa), indicating that the flares occur near the magnetic pole of the 
neutron star. Although the LMC X-4 flares have been observed for many years with various X-ray 
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satellites, no convincing mechanism has been proposed to explain the origin of the flares (e.g., 
Levine et al. 2000). 

On the other hand, LMC X-4 flares have been known as super-Eddington phenomena: some- 
times the flare X-ray luminosity reaches up to '^10^^ ergs s^^ (while its normal state luminosity is 
~2 X 10^^ ergs s~^, comparable to the Eddington luminosity of the ~1.4 Mq neutron star of LMC 
X-4). Because the key physical parameters of LMC X-4 - such as mass, distance, and magnetic field 
- are relatively well determined to be 1.4 ± 0.3 Mq (Levine et al. 1991), 50 ± 10 kpc (e.g., Kovacs 
2000), and ~10^^ G (La Barbera et al. 2001), respectively, the LMC X-4 flares may provide a rare 
example for studying super-Eddington phenomena thoroughly without assuming the key parame- 
ters. This is not the case for most super-Eddington sources owing to relatively large uncertainties in 
their distances and masses. Considering the recent important controversy over the existence of the 
so-called "intermediate-mass black holes" (e.g.. King et al. 2001; Begelman 2001, 2002), for which 
mass (~10^-10^ Mq) is determined under the assumption that the observed X-ray luminosity is 
the Eddington luminosity, the investigation of the origin of the super-Eddington radiation from the 
LMC X-4 flares can also be useful to investigate the origin of the claimed intermediate-mass black 
holes. 

In this paper, we present by far the most extensive observations of the LMC X-4 flares made 
with the Rossi X-ray Timing Explorer (RXTE). We describe the RXTE observations and our data 
analyses in §2; we present the results of timing and spectral analyses in §3. We discuss the possible 
origin of the LMC X-4 flares, as well as their super-Eddington radiation, in §4, We summarize our 
results in §5. 

2. Observations and Data Analyses 

RXTE (Bradt, Rothschild, k Swank 1993) observed two large flaring episodes of LMC X-4 on 
1996 August 19 and 1999 December 19 during its long-exposure (49 and 42 h) observations of the 
source. The HEAsoft^ package (version 5.1) was used for analyzing the data from the Proportional 

Counter Array (PC A; Jahoda et al. 1996). The data within 30 minutes after passages through 
the South Atlantic Anomaly and/or with a higher (> 0.1) electron ratio were ignored. Only the 
data obtained with three or more active Proportional Counter Units (PCUs) were used for timing 
analyses; only the data obtained with all five PCUs were used for spectral analyses. The photon 
arrival times from the Good Xenon mode were transformed to the solar system barycenter using 
the JPL DE400 cphcmcris. The Very Large Event Models of Epoch 3 and 4 were used to subtract 
backgrounds for the data obtained in 1996 and 1999, respectively. The Standard 2 data obtained 
only from the top xenon layers of the PCUs 0, 1, 2, and 4 were used in spectral fits due to 
responsivity problems of PCU 3 (R. Remillard 2002, private communication). Only the spectrum 



^http: / /heasarc.gsfc. nasa.gov/docs / software/lheasoft / 



-4- 



between 2.5-25 keV range was considered for the same reason, and a systematic uncertainty of 1 
% was assumed in spectral fits. A total of 45 data segments (24 from the 1996 observations and 21 
from the 1999 ones) of '^l-h length was obtained, and flares were observed in 4 and 6 segments of 
the two observations. 



3. Results 

3.1. Light curves and Rapid Aperiodic Variability 
3.1.1. Light curves 

Because no reliable ephemeris of the ~30.5-d super-orbital period of LMC X-4 around the two 
observations reported here was available, we used the RXTE All Sky Monitor (ASM; Bradt et al. 
1993) data to examine the phase of the super-orbital motion. Figure 1 shows the ASM light curves 
around the two observations, indicating that they are likely for the high state of the source. 

Figure 2 presents background-subtracted light curves of the PCA data - (a) for the 1996 
observations; (b) for the 1999 ones. The Modified Julian Dates (MJDs) of the start of the two light 
curves are 50314.07905 and 51531.48489, respectively. The binary orbital phase in the upper x-axis 
was determined by the LMC X-4 ephemeris of Levine et al. (2000). The flares are easily identified in 
the first four data segments of Figure 2a and six data segments of Figure 2b as significant increases 
in the photon count rates. Whereas the first half of the data of the 1996 observations shows strong 
flaring activity in the orbital phases (f) ~ 0.26-0.6, the second half shows no noticeable activity at 
the same orbital phases (although some phases arc missing due to data gaps). The LMC X-4 eclipse 
appears around (j) = 1.0 as an apparent decrease in the photon count rates. The 1999 observations 
show flaring activity at ~ 0.4-0.8. The mean photon count rates of the normal state is ~44 and 
^^54 counts/sec/PCU for the 1996 and the 1999 observations, respectively, and that of the eclipse 
is ~13 counts/sec/PCU. (We refer to the normal state in this paper as a state absent of flaring 
activity and eclipse.) 

Figure 3 presents the magnified light curves of the ten data segments containing flares with 
4-s resolution. We shall call them "FLl" to "FLIO" in time order. Figure 3 shows the complex 
flaring activities of LMC X-4 summarized as follows. First, the light curves with strong flares (i.e., 
FLl-5) appear to consist of multiple peaks embedded in the relatively low-intensity regions that 
are still brighter than the normal state; those with weak flares (i.e., FL6-10) show simple, small 
peaks connected through the regions of which brightness is comparable to that of the normal state. 
Secondly, FLl and FL2 clearly show recurrences of intense flares; FL3 and FL5 show a strong 
concentration of intense flares with a time scale of ~2000 s. Finally, while most of the flare peaks 
show symmetric, Gaussian-like profiles, the largest ones (e.g., those at t ~ 3240 s of FL3, t ~ 1770 
s of FL4, and t ~ 2010 s of FL5) appear to show rather asymmetric profiles with a slow decay. 

For quantitative analyses, we performed multi-component Gaussian fits to the flare peaks 
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in Figure 3. Table 1 summarizes parameters of the 68 Gaussian-fitted flare peaks of which peak 
intensity is greater than 90 counts/sec/PCU. The mean peak intensity and full-width-half-maximum 
(FWHM) are 172 ± 103 counts/sec/PCU and 68 ± 31 s, respectively. The median and mean values 
of the reduced chi-square (= Xu) of ^^s are 1.6 and 3.8 ± 5.7, respectively. This indicates that 
while most of the flare peaks are fitted reasonably well (i.e., xl < 2) with a Gaussian function, 
some significantly deviate from a Gaussian. If we exclude the six fiare peaks with Xv greater than 
10 (including the three aforementioned largest fiare peaks that apparently show slow decays), the 
mean Xu decreases to 2.2 ib 1.9. The mean fluence of the flare peaks is (3.1 it 2.9) x 10^*^ ergs, 
with the maximum at 1.9 x 10*^^ ergs. The fluence is the integrated luminosity over the FWHM, 
and the luminosity was computed by spectral fits (see §3.2). (All the errors quoted here represent 
l-cr rms deviation.) 

In Figure 3, the light curves FL2, FL3, and FL4 begin with a normal state, while other light 
curves begin with a flaring state. The mean photon count rates of the first 500 s of FL2, FL3, and 
FL4 are ~35, ^39, and ^40 counts/sec/PCU, respectively, similar to those of the normal state. 
This off'ers an opportunity to examine the light curve transition from a normal states to a flaring 
state. In fact, the magnified views around the beginnings of the flares in FL2, FL3, and FL4 (Figure 
4) identify the existence of the rather flat precursors of the flares lasting for a few minutes just 
before abrupt increases in photon count rates. 

In contrast with the complicated differences in the flare light curves (Figure 3), the pulse 
proflles of FLl-10 show a remarkable similarity (Figure 5), all having simple, sinusoidal profiles. 
Figure 6 compares the pulsed fraction with the peak intensity of the 10 pulse profiles in Figure 5. 
The pulsed fraction is proportional to the peak intensity of the pulse profile with a linear correlation 
coefficient of ~0.76. 

3.1.2. Rapid Aperiodic Variability 

In order to examine the possible existence of any rapid variability associated with the flares, 
we present 45-s light curves of FLl-10 with 0.0625-s resolution around the brightest peak in each 
light curve (Figure 7). The light curves are modulated by ~13.5-s pulsations of the neutron star of 
LMC X-4, and exhibit strong narrow peaks of which intensities are significantly greater than those 
expected from the Poisson noise of the given light curve. The Leahy-normalized power-density 
spectra (PDSs) of the light curves of the central 10-s segments reveal the existence of significant 
powers up to a few hertz, especially for the bright flares such as FL2-5. (Figure 8; We only used 
the central 10-s segments in order to avoid the contamination by the source's pulsations in PDSs. 
Note that the Leahy power owing to a Poisson noise is 2 [Leahy et al. 1983]).) We investigated 
this possible correlation between the existence of the signiflcant powers in the PDSs and the flare 
intensity more as follows. First, we divided the FLl-10 light curves (Figure 3) into segments of 
which lengths are equal to the pulsational period (~13.5 s) of LMC X-4 determined by Fourier 
transformations of the light curves. Next, we computed the Leahy-normalized PDSs of the all 



-6- 



~13.5-s segments, as well as their integrated intensities. The correlation between the flare intensity 

and the existence of the significant powers in PDSs is evident in Figure 9, where we compare the 
distributions of the integrated intensities with those of the significant powers (> 10) (in PDSs) of 
the ~13.5-s segments. This confirms that the rapid aperiodic variability during the flares increases 
along with the flare intensity, which sometimes appears as strong narrow peaks in the light curves 
(Figure 7). 

3.2. Spectral Evolution of the Flares 

3.2.1. Softness Ratio Distribution 

We investigated the spectral evolution of the flares via the softness ratio, which we defined 
to be the ratio of the soft X-ray (2-8 keV) photon count rates to those of the hard X-ray (10-20 
keV). Figure 10 compares the total photon count rates (integrated over the 2-25 keV range) of the 
flares with the softness ratios obtained with 32-s resolution, identifying a strong linear correlation 
between them with a linear correlation coefficient of ~0.95. 

However, Figure 11a, which is the same as Figure 10 but obtained with 0.0625-s resolution, 
shows that the strong correlation identified with 32-s resolution disappears when time resolution 
improves - the linear correlation coefficient decreases to ~0.30. The data points with large photon 
count rates (so with small Poisson noise) show more significant deviations from the linear correlation 
(identified with 32-s resolution) than those with small photon count rates, indicating that the 
deviation is not simply owing to the increased Poisson noise caused by the improved time resolution. 
We investigated this further via comparison of the correlations obtained from two separate groups 
of the data: one is the "faint group" having small total photon count rates; the other is the "bright 
group" having large total photon count rates. We chose the faint group as the group of data points 
of which total photon count rates belong to the lower 90 %; the bright group of which total photon 
count rates belong to the upper 1 %. The dotted and dashed lines in Figure 11a represent the 
correlations obtained in the faint and bright groups, respectively. The correlations obtained in the 
two groups are quite distinctive: while the faint group (dotted line) shows a positive correlation 
with a slope similar to that obtained with 32-s resolution (Figure 10), the bright group (dashed 
line) docs not show any apparent correlation, indicating that the strong linear correlation between 
the fiare intensities and the softness ratios obtained with 32-s resolution does not hold for this 
component (i.e., the narrow, spiky peaks in Figure 7). 

In order to examine the effect of the increased Poisson noise caused by the improved time 
resolution in Figure 11a more thoroughly, we used a Monte Carlo simulation by re-sampling the 
each data point obtained with 32-s resolution (Figure 10) with 0.0625-s resolution. We re-sampled 
the 32-s resolution data points using a Gaussian distribution with the la standard deviation equal 
to the Poisson noise of 0.0625-s resolution. We performed the simulation 1000 times, and present 
a typical result in Figure lib. Although the correlation coefficients of the three lines representing 
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the all data points (solid line, ~0.20), the faint group (dotted line, ~0.09), and the bright group 
(dashed line, ~0.25) are much smaller than that obtained with 32-s resolution (~0.95), they show 
very similar slopes to that obtained with 32-s resolution and no apparent difference among them. 
This contrasts with the real data (Figure 11a), demonstrating that the simulated data are unable 
to produce the decreases in the softness ratios of the bright group of the real data (while the 
correlation coefficients of the simulated data decrease due to the increased Poisson noise caused by 
the improved time resolution, 0.0625 s). Therefore, the absence of the linear correlation between 
the softness ratios and the total photon count rates in the bright group of the real data seems to 
be of statistical significance. Table 2 summarizes the correlation coefficients and slopes of the lines 
in Figures 10 and 11. 

The systematic difference in the softness ratios of the bright and faint groups of the real data 
is more clearly identified in Figure 12, where we compare the binned softness ratios of the real data 
(dotted line) in Figure 11a with those of the simulated data (solid line) in Figure lib. For this, we 
made 100 bins of the total photon count rates (each bin with ~1230 data points), and computed 
the mean softness ratios in each bin. As we see in Figure 12, the softness ratios of the real data 
with small total photon count rates (< 300 counts/sec/PCU) are almost identical to those of the 
simulated data; however, the real data show a clear systematic decrease in the softness ratios as 
the total photon count rates increase (> 300 counts/sec/PCU). About 85 % of the data have total 
photon count rates less than 300 counts/sec/PCU. 

3.2.2. Spectral Fits 

In order to study the spectral evolution of the flares more thoroughly, we obtained flare spec- 
tra over 32-s interval, and fitted them to a model spectrum consisting of a power law, a black 
body, and a Gaussian component, as well as a component for photoelectric absorption due to the 
intervening interstellar matter. The power law and Gaussian components represent non-thermal 
magnetospheric emission and iron line emission, respectively, which is typical of X-ray pulsar (e.g., 
White, Swank, & Holt 1983). The black body component represents the thermal emission asso- 
ciated with the neutron star polar cap and the accretion columns during the flares. The model 
spectrum is given by 

I{E) = e^p[-a{E)Nti] x {IpE"" + BB{E) + G{E)} photons cm^^ s"^ keV"^ (1) 

where E is the photon energy in kcV; o-(E) and A^h are the photoelectric cross section and the 
hydrogen nuclei column density of the intervening matter; Ip and a are the norm and the power- 
law index of the power-law component. BB{E) is the black body function of temperature kT; 
G{E) is the Gaussian function for the iron line emission with its center at Eq. The lower limit of 
A^H was fixed to be 3 x 10^° cm"^ based on the previous results (Woo et al. 1995, and references 
therein). We obtained the mean of 0.99 ± 0.22 with ~43 degrees of freedom from spectral fits 
between 2.5 and 25 keV. 
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Figure 13 presents the distributions of the resulting power index (a, 13a) and the black body 
temperature (kT, 13b) along with the flux obtained from the model fits in the 2.5-25 keV range. 
The power index increases as the flux increases; the black body temperature decreases on the 
contrary. The possible saturation in the decrease of the black body temperature may be due to 
the lower energy limit of the fits (= 2.5 keV). Figure 13 shows that both the power-law and black 
body components soften as the flare intensity increases, which is consistent with the softness ratio 
distribution obtained with 32-s resolution (Figure 10). (We tried many different model spectra, 
e.g., adding a high-energy cutoff to the power-law component, subtracting the Gaussian component 
or the black body component, and fixing the hydrogen nuclei column density. The fits to these trial 
spectra gave us a moderately reasonable xt (< 2)- The trends identified in Figure 13 always existed 
in the spectral fits to these trial model spectra. The fits to the model spectrum of Eqn. [1] yielded 
the least xt-) Table 3 lists the luminosities of the precursors (Figure 4) and the brightest flare peak 
of each FLl-10 (Figure 3), as well as the fluence of FLl-10. We used a distance of 50 kpc to LMC 
X-4 to estimate the luminosity. Figure 14 presents a comparison between the best-fit model spectra 
(solid lines) and the observed ones (crosses) of the brightest fiare peak of each FLl-10. Table 4 
summarizes the best-fit parameters of Figure 14. 

We compared the spectrum from the normal state with that of the brightest flare peak in all 
FLl-10 (i.e., the peak at t ~ 1770 s of FL4 in Figure 3). We used three different model spectra: 
(1) the model spectrum given by the Eqn. (1), (2) a model spectrum consisting only of a power-law 
component and a Gaussian component, (3) the same model spectrum of (2) but with a high-energy 
cutoff for the power-law component. For the case (3), the model spectrum is given by 

I{E) = e^p[-a{E)Nn] x {IpE''' fhi{E) + G{E)} photons cm-^ s"^ keV"^ (2) 

where 

J 1, E <Ec 

/'^^(^) = |exp(-(^), E>E,, 

represents the high-energy cutoff for the power-law component. Figure 15 compares the three best- 
fit model spectra (solid lines) with the observed ones (crosses), and Table 5 summarizes the best-fit 
parameters. All the spectral fits to the three model spectra gave us acceptable xt (< 2). We noted 
that the fitted values for the high-energy cutoff component of the normal state were very close to 
those previously obtained (16.1 it 0.4 and 35.6 it 6.7 keV; Woo et al. 1996), whereas those for 
the flare peak were considerably different. In order to examine the significance of this change in 
the high-energy cutoff component, we fit the flare peak spectrum to the model spectrum given by 
the Eqn. (2) with the high-energy cutoff component flxed to the values obtained from the normal 
state: Ec = 16.1 keV and E{ = 27 A kcV. The resulting best-fit parameters arc: A'h = (6.5 it 0.7) 
X 10^2 cm~2^ a = 3.3 ± 0.2, and Eq = ^.S^^qI keV (with xl - 1-3), similar to the values in Table 
5. Therefore, the changes in the high-energy cutoff component between the normal state and flares 
do not seem to be significant. The mean luminosity of the normal state luminosity is ~1. 5-2.0 x 
10^^ ergs s~^ 
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4. Discussion 

4.1. Possible Origin of tlie Flares 

4.1.1. Thermonuclear Burning 

Levine et al. (2000) discussed the possibility of carbon burning under the neutron star surface 
as the origin of the LMC X-4 flares. According to their discussion, the nuclear burning scenario can 
roughly explain several features of the flares, including the luminosity, the daily recurrence time, 
the rapid rise and gradual decay, and the spectral softening. Shortcomings are an absence of any 
mechanism which can explain the '^100-200 s substructure and the lack of enough carbon material 
needed for powering the flares when the accreting material is mostly H-rich material (Brown & 
Bildsten 1998). Although it might be possible to attribute the ~100-200 s substructure to an 
accretion disk origin (Moon &; Eikenberry 2001), the results of this paper are not supportive of the 
nuclear burning mechanism as the origin of the flares. First, the precursors found in this paper 
(Figure 4) do not reconcile readily with the nuclear burning scenario, because the higher accretion 
rate during the precursor on the polar cap region would suppress the thermonuclear instability (e.g., 
Bildsten & Brown 1997; Brown &; Bildsten 1998). We note that Type I X-ray bursts, which are 
known to be caused by nuclear flashes on the neutron star surface in low-mass X-ray binaries, have 
occasionally been observed with precursors of a few seconds (e.g., Lewin, Vacca, & Basinska 1984). 
However, the burst precursors are distinctive from the precursors of the LMC X-4 flares of this 
paper. The burst precursor is generally followed by a normal state of a few seconds due to a rapid 
expansion and subsequent cooling of the neutron star photosphere, which is clearly different from 
the LMC X-4 flares. In addition, while many Type I X-ray bursts have been observed with spectral 
softening during the decay, the LMC X-4 flares show the opposite phenomena: the spectrum is 
hardening during the decay. The very strong magnetic field strength of LMC X-4 (^10^^ G; La 
Barber a et al. 2001) also maJies the nuclear burning scenario less likely - no X-ray pulsars have 
been observed with Type I X-ray bursts except for rare accreting millisecond pulsars (e.g., SAX 
J1808. 4-3658; in 't Zand et al. 1998). Finally, the observed spectra of the flares are far from 
thermal, although the thermal spectra from a highly magnetized object may be very different from 
those of non-magnetized objects. Therefore, it is very difficult to explain the LMC X-4 flares in 
terms of the nuclear burning scenario. 



4.1.2. Viscous Instability 

We note that X-ray bursts/flares from the X-ray pulsars GRO J1744-28 and SMC X-1 have 
been explained as results of viscous instabilities in the accretion disk (Cannizzo 1996, 1997; Li &i 
van den Hcuvel 1997; Moon et al. 2002). If the radiation pressure is comparable to the gas pressure 
around the inner-disk radius, the accretion disk is susceptible to a viscous instability. Under the 
classical a-disk configuration (Shakura & Sunyaev 1973), this requires that the transition layer 
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between the radiation-pressure-dominated and the gas-pressure-dominated regions be comparable 
to the inner-disk radius. In order to satisfy this condition, the magnetic moment of the neutron 
star is expected to be (Shapiro & Teukolsky 1983) 

M30 - 4 X 10-=^ MiV'/^ (3) 

where /Ltao is the magnetic moment of a neutron star in unit of lO^'^ G cm^, and M17 is the mass 
accretion rate in unit of 10^^ g s~^. We assumed 1.4 Mq and 10 km for the mass and radius of 
the neutron star, respectively, and obtained /X30 — 0.3 with M17 2± 10. This is significantly smaller 
than the value estimated from observations (/X30 c± 11; La Barbera et al. 2001), indicating that 
the viscous instability is difficult to grow in the accretion disk around LMC X-4 with the a-disk 
configuration. This also confirms that it is unlikely that radiation-pressure-dominated regions can 
exist around highly-magnetized neutron stars (i.e., X-ray pulsars). 



4.1.3. Plasma Rayleigh- Taylor Instability and Stellar Winds 

LMC X-4 is known to be in a spin equilibrium state (Naranan et al. 1985; Woo et al. 1995; 
La Barbera 2001) in which the magnetospheric radius is comparable to the corotation radius of 
the source. In this case, it has been suggested that the plasma Rayleigh-Taylor instability (which 
is known as the Kruskal-Schwarzschild instability) can develop to result in X-ray bursts from the 
neutron star X-ray binaries (Baan 1977, 1979). According to this scenario, accreting matter can be 
accumulated outside the magnctosphere when the neutron star spins faster than the inner accretion 
disk owing to the centrifugal barrier of the magnctosphere. Once the density of the accumulating 
matter exceeds the critical value, the effective gravity becomes larger, driving the magnctosphere 
inside the corotation radius, and an abrupt increase in the accretion onto the neutron star can occur. 
Based on this mechanism, Apparao (1991) explained the origin of the six flares observed from EXO 
2030+375. We note that the flares of EXO 2030+375 share similar properties with the LMC X-4 
flares, including the ~3-4 h flare recurrence time, the ~l-2 h flare duration, the high fluences of 
the flares (> lO^'^ ergs), the very similar precursors, the rapid rise with gradual decay, and the 
absence of any significant spectral change (Parmar et al. 1989; Parmar, White, Sz Stella 1989). 
Similar to LMC X-4, ~10^^ G magnetic field strength is expected from EXO 2030+375 (Apparao 
1991). These similar properties are suggestive that the same mechanism may be responsible for the 
flares from the both objects. 

Finally, we investigate the possibility that the LMC X-4 flares are caused by non-uniform 
stellar winds from the 0-type companion. One noticeable result is that, based on the estimation 
of stellar mass loss rate, Boroson et al. (1999) found that the gravitational capture of stellar 
winds could power ~3 x 10^'' ergs s~^ luminosity of LMC X-4, which is comparable to the excess 
luminosities of the flare precursors over the normal states (Table 3). In fact, several early type 
stars have been observed with non-uniform stellar winds of ~2-3 h periods (Vogt &; Penrod 1983; 
Kudritzki Sz Puis 2000) , which is similar to the time interval of the three consecutive flares of the 
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LMC X-4 (FLl-3; Figure 2 and 3). (The flares FL6-10 are also consecutive, but it is difficult to 
estimate flare recurrence time from these flares.) Additionally, the roughly estimated recurrence 
period of the flaring episode of LMC X-4 is compatible to its orbital period, '^1.4 d. Based on these, 
one can imagine a very speculative scenario in which the LMC X-4 flare precursors are caused by 
the capture of material in the non-uniform stellar winds from the companion star related to binary 
orbital motion, and then the precursor induces instabilities in the accretion disk (such as the plasma 
Rayleigh- Taylor instability) by satisfying the conditions for the instability to develop. 

4.2. Softness Ratio Distributions and Inhomogeneous Polar Cap During 

Super-Eddington Flares 

4-2.1. Softness Ratio Distributions and Inhomogeneous Polar Cap Accretion Flows 

As identified by the analysis with 0.0625-s resolution, the LMC X-4 flares appear to consist of 
two distinctive components. The first one, which is the major component corresponding to ~85 % 
of the total data, has relatively small total photon count rates (< 300 counts/sec/PCU); the second 
one, which is the minor component corresponding to the remaining ~15 %, has large total photon 
count rates (> 300 counts/scc/PCU). The softness ratio distribution of the major component is 
almost equal to that obtained with 32-s resolution, showing a strong linear correlation between the 
softness ratio and flare intensity (Figure 10 and 12). The minor component also shows an overall 
proportionality between them; however, it has a distinctive feature showing a much slower increase 
in the softness ratio, as the flare intensity increases, than the major component (Figure 12). The 
minor component often appears as very narrow, spiky peaks in the light curves (Figure 7), and is 
also responsible for the increased rapid aperiodic variability in PDSs (Figure 8). We shall call the 
major component the "broad, normal" component; the minor component the "narrow, spiky" one. 
With 32-s resolution, the narrow, spiky component is dominated by the broad, normal component, 
so that it is impossible to identify them. 

Independent of the time resolution used, the LMC X-4 flare spectra become softer as the flare 
intensities increase (Figure 10, 12, and 13). For the case of high-luminosity sources (> 10^^ ergs 
for example), the accretion column is expected to be optically thick to its own X-ray radiation and 
to re-radiate in soft X-rays (Basko & Sunyaev 1976). Because the number of reprocessings (i.e., 
scatterings/absorptions) of photons are expected to be proportional to the flare intensities (i.c, the 
accretion rates), it seems natural to attribute the origin of the broad, normal component to the 
increased optical depths of the accretion columns during the flares, which can explain the spectral 
softening along with the flare intensities. (See Basko & Sunyaev [1976] for detailed mechanisms 
responsible for the spectral softening. We assume here that the LMC X-4 flares are caused by the 
abruptly increased accretion rates possibly caused by an accretion disk instability [see §4.1].) 

This mechanism can also account for the proportionality between the softness ratios and flare 
intensities of the narrow, spiky component. However, it is not readily capable of explaining the 
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much slower increase in the softness ratios as the flare intensities increase (Figure 12), indicating 
that a different mechanism is responsible for it. Wc propose here that density inhomogeneity in 
the accretion columns onto the neutron star polar cap during the flares may be responsible for the 
observed softness ratio distributions. If the accretion column is inhomogeneous, it is possible that 
some photons escape it without suffering as many reprocessings as expected by the flare intensities. 
As indicated by the increase of the rapid aperiodic variability (Figure 9), the inhomogeneity in 
the accretion columns is most likely proportional to the flare intensities, suggesting that the prob- 
ability of the leakage of photons is higher for the intcnscr flares. In this case, the photons from 
more intense flares are expected to have a greater chance of leaking through the inhomogeneous 
accretion columns. Therefore, while the overall softness ratios are still proportional to the flare 
intensities, the intense flares are expected to show a slower increase in the softness ratio as the 
flare intensities increase - consistent with the observed softness ratio distributions (Figure 10 and 
12). In contrasts with the large scatter in the softness ratios of photons having similar total photon 
count rates (Figure 10), the statistically averaged softness ratios show a well-defined systematic 
behavior (Figure 12). This is also consistent with the argument given above, because it is based on 
the probability of the leakage of photons, which can best appear as a statistical behaviour. In con- 
clusion, the narrow, spiky peaks in the light curves (Figure 7) seem to represent the photons leaking 
through the inhomogeneous accretion columns; the rapid aperiodic variability proportional to the 
flare intensity (Figures 8 and 9) is most likely related to the scale of the density inhomogeneity in 
the accretion columns. 

For bright X-ray pulsars, the previous studies have shown that low-density "photon bubbles" 
form at the settling mound of the accretion column above a polar cap (Arons 1992), because the 
relatively heavy accretion column is supported against gravity by radiation. If the X-ray radiation 
is strong enough (L > 10^^ ergs s~^), a radiative shock forms at the bottom of the accretion column 
(Basko & Sunyaev 1976; Wang & Frank 1981; White et al. 1995). Numerical studies (Klein et 
al. 1996a) have shown a variation in the number of photon bubbles if accretion continues onto 
the neutron star polar cap, resulting in quasi-periodic variations in the luminosity by changing the 
efficiency of radiation transport in the accretion column. These are known as the "photon-bubble 
oscillations" (PBOs). Klein et al. (1996b) have found the relation, z^pbo L~^-^^^, between the 
frequency of the PBOs and luminosity in the numerical simulation for the narrow polar cap model, 
assuming 9c = 0.1, B = 3 x 10^^ G, and L = (0.07-2) x 10^^ ergs s"^ for the polar cap radius, the 
neutron star magnetic field, and the luminosity, respectively. This narrow polar cap model predicts 
fast (> 1000 Hz) PBOs. Due to the increased photon diffusion time, Klein et al. (1996b) have 
also found that the frequency of PBOs decreases if the polar cap size increases. According to their 
results on the larger polar cap model (^c = 0.4), ~85 Hz PBOs are possible for an X-ray pulsar 
with L = 3 X 10^^ ergs and B = 3 x 10^^ G. Because no previous study is directly applicable 
to the extreme conditions of the LMC X-4 flares studied in this paper (L ~ 10^^ erg s~^; B ~ 10^^ 
G), we extended the results of Klein et al. (1996b) to the LMC X-4 flares. If we use the upBO oc 
^-0.617 scaling law obtained from the narrow polar cap model and i'pbo = 85 Hz from the larger 
polar cap model of L = 3 x 10^^ ergs s~^, the expected PBO frequency of the LMC X-4 flares is 
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~10 Hz for L ~ 10 ergs s~ , comparable to the observed frequencies of the aperiodic variabilities 
during the flares (Figure 8). We need a more through numerical simulation directly applicable to 
the LMC X-4 flares to investigate the possible relations between the rapid aperiodic variability 
discovered in this paper and PBOs. 

4-2.2. Origin of the Super-Eddington Radiation of the Flares 

The flare luminosities obtained in this paper are often well in excess of the Eddington luminosity 
of the 1.4 Mq neutron star in LMC X-4, ~2 x 10^^ ergs s"-*^. However, the Eddington luminosity 
is the upper limit that a steady, isotropic, homogeneous, and fully ionized atmosphere can support, 
which is not the case for LMC X-4 primarily because that most emission is through the neutron 
star polar cap. The size of the neutron star polar cap is expected to be much smaller than the 
total neutron star surface (i.e., < 10 %). Therefore, if all the observed emission is through the 
polar cap, this indicates that sometimes LMC X-4 flare luminosities are above the Eddington limit 
by more than two orders of magnitude. However, because a significant fraction of the emission 
may be transported as fan beams through the sides of accretion columns above the polar cap 
(due to the increased optical depths along the direction of the magnetic field), it is not clear if 
the polar cap emission can really decrease the Eddington limit by such magnitude. On the other 
hand, the photon scattering cross-section in a strongly magnetized plasma has been known to be 
considerably reduced along the magnetic field lines for the photons of which energy is much lower 
than the cyclotron energy (White, Nagese, & Parmar 1995, and references therein). According to 
the results of Paczyhski (1992), the strong magnetic field (~10^^ G) of LMC X-4 can render its X-ray 
luminosity up to ~50 times of the Eddington limit at the neutron star sTirfacc by reducing Thomson 
and Compton cross sections. However, this effect is expected to be substantially diminished during 
the LMC X-4 flares. This is because, during the flares, radiative shocks can easily form at the 
bottom of the accretion columns, supporting the accretion columns above the neutron star polar 
cap (Basko & Sunyaev 1976; Wang &; Prank 1981; White et al. 1995), which would significantly 
decrease the effect of the magnetic field on the photon cross sections (note that the magnetic field 
is decaying 1/r^). Therefore, we regard the LMC X-4 flares as super-Eddington phenomena. 

What then is the origin of the super-Eddington radiation of the LMC X-4 flares ? We suggest 
that the inhomogeneous accretion columns during the flares, which were already proposed to explain 
the softness ratio distributions and rapid aperiodic variability (§4.2.1), are also responsible for the 
super-Eddington radiation. Recently, some authors have proposed that a magnetized, radiation 
pressure-dominated atmosphere (and an accretion disk) is susceptible to a density inhomogeneity 
(e.g., Shaviv 1998, 2000; Begelman 2001, 2002), most likely caused by the non-linear development 
of a "photon-bubble instability" (Arons 1992; Gammie 1998). In this mechanism, once the matter 
is trapped in the magnetic field by radiation pressure, the magnetic field prevents the trapped 
matter from spreading sideways, maintaining the density inhomogeneity. The net flux through 
this inhomogeneous atmosphere can easily exceed the Eddington limit by a large factor primarily 
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through the low-density (= low-opacity) regions (Shavie 1998, 2000; Begelman 2001). We note 
that the LMC X-4 flares satisfy the conditions for the development of the density inhomogeneity 
in the accretion columns easily owing to their high luminosities and the strong magnetic field of 
the source. In addition, because a radiative shock can readily develop above the polar cap due to 
the high luminosities of the LMC X-4 flares, the "photon bubbles" are expected to form at the 
bottom of the accretion columns during the flares, which can be a seed for the proposed density 
inhomogeneity. 

4.3. Pulse Profiles 

During the LMC X-4 flares, the accretion columns are expected to be highly opaque. In this 
case, one might expect that a significant fraction of the flare emission is transported via the fan 
beams in the direction perpendicular to the magnetic field due to the relatively smaller optical 
depths in this direction than the direction parallel to the magnetic field. The increased fan beams 
then can result in substantial gravitational bendings of the X-ray photons from the neutron star 
polar cap (e.g., Brainerd k, Meszaros 1991), which may be responsible for the broad pulse profiles 
of the LMC X-4 flares. On the other hand, the simple sinusoidal, in-phase pulse profiles may 
indicate that pulse profiles are simply dominated by one mechanism, i.e., the flares, probably at 
one magnetic pole. 

5. Summary and Conclusion 

In this paper, we have analyzed two extensive observations of the LMC X-4 flares made with 
RXTE. We summarize our main results as follows. (1) The flare light curves reveal the existence of 
the fiat precursors preceding the flares, with the complex activity of the LMC X-4 flares. (2) The 
FWHM and mean fluence of the flares are 68 ± 31 s and (3.1 ± 2.9) x 10^° ergs, respectively. (3) 
The pulse profiles during the flares are broad, sinusoidal, and the pulsed fractions increase with the 
flare intensities. (4) With 0.0625-s resolution, narrow, spiky peaks appear in the light curves, so 
does the rapid aperiodic variability in PDSs. Both intensify along with the flare intensities. (6) The 
softness ratios and spectral fits show that the flare spectra soften as the intensities increase. (7) 
With 0.0625-s resolution, while the photons having count rates less than 300 counts/sec/PCU show 
the softness ratio distributions almost equal to those obtained with 32-s resolution, the photons 
having larger photon count rates show a much slower increase in the softness ratios as the flare 
intensities increase. (8) The peak luminosity of the flares reaches up to ~2.1 x 10^^ ergs s~^, well 
above the Eddington limit for an 1.4 Mq neutron star. 

We find it difficult to explain the observed properties of the LMC X-4 flares in terms of 
the thermonuclear instability on the neutron star surface or the viscous instability at the inner 
accretion disk. We consider the possibility that the plasma Rayleigh- Taylor instability in the 
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accretion disk may be responsible for the flares. We propose that increased optical depths and 
density inhomogencity in the accretion columns onto the neutron star polar cap during the LMC 
X-4 flares arc responsible for the observed properties, including the super-Eddington luminosities of 
the flares. If confirmed, this is the first observational evidence for the inhomogeneous atmosphere 
responsible for super-Eddington radiation. The broad pulse profiles of the LMC X-4 flares may be 
caused by substantial gravitational bending of the flare emission transported via fan beams. 
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Table 1. Parameters of the LMC X-4 Flare Peaks 
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Table 1 — Continued 
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Table 1 — Continued 



Flare 


position 


peak* 


FWHM 


dof 




fluence'^ 




(sec) 




(sec) 






(10^9 ergs) 


FLIO 


2808 


103 


51 


7 


1.5 


29 


FLIO 


3261 


136 


58 


7 


1.9 


37 


FLIO 


3600 


135 


67 


9 


1.7 


39 


FLIO 


4009 


155 


54 


7 


1.1 


23 


FLIO 


4461 


147 


51 


6 


4.6 


24 


FLIO 


4565 


95 


78 


10 


1.2 


23 



*In unit of counts/sec/PCU. 
^Caculated at 2.5-25 keV. 
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Table 2. Correlation Coefficients and Slopes of the Lines in Figures 10 and 11 



Resolution 


All data 


Faint group 


Bright group 


32 s^ 


0.95, 1.4 X 10-2 (240) 






0.0625 


0.30, 7.8 X 10-3 (122880) 


0.19, 1.1 X 10-2 (110592) 


-0.01, -2.4 X lO-"* (1228) 


0.0625 s^ 


0.20(13), 1.0(0) X 10-2 


0.09(6), 9.8(0) X 10-3 


0.25(4), 8.2(6) X 10-^ 



^From the real data. The numbers in parentheses are the number of data points. 

'^From the Monte Carlo simulation. The numbers in parentheses correspond to la deviations at 
the last quoted digit. 
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Table 3. Luminosity and Fluence of FLl-10 



Name precursor peak fluence^ 

(10^^ ergs s-^) (10^^ ergs s"^) (10^^ ergs) 



FLl 




6.1 




1.0 


0.7 


± 


0.2 


FL2 


1.8 ± 0.5 


13.0 




1.6 


1.4 




0.3 


FL3 


1.9 ± 0.4 


12.9 




1.7 


1.1 




0.3 


FL4 


2.0 ± 0.3 


21.4 




2.6 


1.9 




0.4 


FL5 




4.8 




0.2 


0.8 




0.3 


FL6 




2.3 




0.3 


0.5 


± 


0.2 


FL7 




1.9 




0.3 


0.5 




0.2 


FL8 




1.7 




0.4 


0.5 


± 


0.3 


FL9 




3.9 




0.8 


1.0 




0.3 


FLIO 




G.l 


± 


0.9 


1.1 




0.3 



^Including the normal state. 

Note. — The errors correspond to 1 a uncertanties, 
and the energy range is 2.5-25 keV 
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Table 4. Best-fit Parameters of the Brightest Peaks in FLl-10 



Name A^h 


a 


kT flux^ xl 


(10^2 cm-2) 




(keV) (10-'' ergs s"^) 



FLl 


^•^+1.2 


1, 


O-0.2 
■°+0.2 


2.3: 


FL2 


0(5-1.0 


2, 


c:-0.2 
•^+0.1 


1.7: 


FL3 


9 c-0.9 


2, 


•^+0.1 


1.6: 


FL4 


4 2"°-^ 


3 


Q-0.2 
■'^+0.1 


1.4: 


FL5 




2, 


r;-0.1 
•^+0.1 


1.5: 


FL6 


1.8^ 


1, 




2.3: 


FL7 


1.8^ 


2, 


0-0.5 
■^+0.4 


4.2: 


FL8 


-■-•^+1.1 


2, 


t;-0.5 
•^+0.5 


4.7: 


FL9 


1.8^ 


2, 


c:-0.6 
■^+0.4 


5.3: 


FLIO 


1 1-1° 
-'-•-'-+1.1 


1, 


Q-0.3 
■^+0.4 


3.9: 



-0.3 
+0.5 
-0.1 
+0.1 
-0.1 
+0.2 
-0.1 
+0.1 
-0.1 
+0.1 
-0.2 
+1.3 
-0.4 
+0.3 
-0.4 
+0.5 
-0.5 
+0.5 
-0.6 
+0.5 



2.0 




0.3 


0.9 


4.3 




0.5 


1.2 


4.2 




0.6 


1.4 


7.2 




0.9 


1.1 


1.6 




0.1 


1.6 


0.8 




0.1 


1.4 


0.6 




0.1 


1.2 


0.6 




0.1 


0.8 


1.6 




0.3 


1.3 


2.0 




0.3 


1.4 



'^In the 2.5-25 keV energy range. 
^The values are fixed in the fits. 



Note. — 



The errors correspond to 1 a uncertanties. 



Table 5. Best-fit Parameters of the Spectral Fits 



Parameters iVn a kT Ef Eq flux^ xl 

(10^2 cm-2) (keV) (keV) (keV) (keV) 



Model 1^ 

Flare Peak 4.8"?;^ 3.1+^;^ 1.4;{]:J 
Normal State 4.5-2:8 ^.l"";^ G.O^^j 



Modol 2' 



Flare Peak 5.7;gj ^-Koi •■■ '■■ '■■ 5- ''+0:3 ^.2 1.4 

Normal State 1.9;^;^ 0.981^^:^^ • • • • • • • • • 7.3;^;^ 0.69 1.7 



Model 3^^ 



Flare Peak 2.9;^;^ 2.0 J;^ ••• 3.6;|| 8.3+Jj 5.9+[]:^ 7.1 1.0 
Normal State 0.60;?;^ O-Sl^oli " " " '^^■'^+2.5 27.4;i|j 6.3^0;^ 0.54 1.5 



'^In unit of of 10 ^ ergs cm ^ s ^. 
•^See the Eqn. (1). 
'^Same as the Eqn. (1), but without the black body component. 
^SeetheEqn. (2). 



6.3-[]J 7.1 1.0 
7.27^:^ 7.1 1.4 




Fig. 1. — The RXTE All Sky Monitor data of LMC X-4 around the two observations reported in 
this paper - (a) for the 1996 August observations and (b) for the 1999 December ones. The arrows 
show the start times of the two observations, and the error bars correspond to la uncertainties. 
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Fig. 2. — Background-subtracted light curves of LMC X-4 in the 2-25 keV range with 32-s resolu- 
tion: (a) for the 1996 observations and (b) for the 1999 ones. Upper x-axes represent orbital phases 
of LMC X-4. 
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Fig. 3. — Same as Figure 2, but only for the data segments containing the flares with 4-s rcsohition. 
(a)-(j) are for FLI-FLIO. The inserted numbers represent the time elapsed from the start of the 
two observations, i.e., the 1996 ones for FLl-4 and the 1999 ones for FL5-10. 
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Fig. 4. — Magnified light curves of FL2, FL3, and FL4 around the beginnings of the flares. The 
inserted numbers represent the time offset from the start of each light curve in Figure 3. 
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Fig. 5. — Pulse profiles of FLl-10 in Figure 3. The phase is adjusted to have the maximum 
photon count rates, and the number of bins is 20. The inserted numbers represent pulsed fractions. 
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Fig. 6. — Comparison between the pulsed fractions and the peak intensities of the pulse profiles in 
Figure 5. 
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Fig. 7. — Same as Figure 3, but only for the 45-s data segments around the brightest peak in 
each hght curve with 0.0625-s resolution. The inserted number represents the time elapsed from 
the start of each light curve; the length of the vertical bars represents the magnitude of the mean 
Poisson noise. 
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Fig. 8. — Leahy-normalized PDSs of the 10-s segments around the brightest peak in FLl-10. 
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Fig. 9. — Distributions of the integrated (over '~13.5-s interval) flare intensities (solid lines; in 
arbitrary logarithmic intensity scale) and signiflcant Leahy powers (> 10; dots). 
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Fig. 10. — Distribution of the softness ratios computed with 32-s resolution. The sohd hne repre- 
sents the linear correlation between the softness ratios and the total photon count rates. The line 
is extended to be compared with Figure 11. 
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Fig. 11. — (a) Same as Figure 10, but with 0.0625-s resoluton. The dotted and dashed hnes 
represent the correlations between the softness ratios and the total photon count rates found in the 
lower 90 % and the upper 1 % data points, respectively. The lines are extended to the full extents 
of the figures to be seen more clearly, (b) Same as Figure (a), but obtained from the Monte Carlo 
simulation. 
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Fig. 12. — Distributions of the softness ratios integrated in 100 total intensity bins. The sohd hne 
is for the real data (Figure 11a); the dotted line is for the simulated data (Figure lib). 
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Fig. 13. — Comparisons of the fitted power index (a) and the black body temperature (b) of the 
flares with the flux at 2.5-25 keV. 
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Fig. 14. — Comparisons between the observed spectra (crosses) and the best-fit spectra (soHd hnes) 
of the flare peaks in FLl-10. 
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Fig. 15. — Same as Figure 14, but for the spectra from the brightest peak in the ah flares (a, c 
and e), and for the spectra typical of the normal state (b, d, and f). (a) and (b) are for model 
spectrum from the Eqn. (1), (c) and (d) are the same with (a) and (b) but without the black body 
component, and (e) and (f) are for the Eqn. (2). 
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